Defect Relaxation and Coarsening Exponents 
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The coarsening exponents describing the growth of long-range order in systems quenched from 
a disordered to an ordered phase are discussed in terms of the decay rate u)(k) for the relaxation 
of a distortion of wavevector k applied to a topological defect. For systems described by order 
parameters with Z(2) ('Ising') and 0(2) ('XY') symmetry the appropriate defects are domain walls 
and vortex lines respectively. From tot ~ k z , we infer L(t) ~ f. 1//z for the coarsening scale, with 
the assumption that defect relaxation provides the dominant coarsening mechanism. The 0(2) case 
requires careful discussion due to infrared divergences associated with the far field of a vortex line. 
Conserved, non-conserved, and 'intermediate' dynamics are considered, with either short-range or 
long-range interactions. In all cases the results agree with an earlier 'energy scaling' analysis. 



Ginzburg-Landau (TDGL) equation (see section II) and 
the Cahn-Hilliard equation (section III). 

Powerful though the energy-scaling method is, the ar- 
guments are not completely straightforward and some 
care is required in its application (see Q] for a detailed 
discussion). In this paper, therefore, we will discuss a 
much 'cleaner' approach that addresses the dynamics of 
the topological defects in a more direct way. Specifically 
we consider a single plane domain wall, or straight vor- 
tex line, and apply a periodic perturbation of wavevector 
k. The perturbation then relaxes away, with asymptotic 
relaxation rate uj(k). A central assumption, as before, is 
the validity of the scaling hypothesis. If there is a sin- 
gle characteristic length scale L(t), then the relaxation 
dynamics of the defect should be reflected in the coars- 
ening dynamics, through L/L ~ u>(k) with k ~ 1/L. So 
if uj(k) ~ k z for k — > 0, we expect L(t) ~ t x l z . More 
simply, one can set u) ~ 1/t and k ~ 1/L to obtain the 
same result. 

This approach has been used before for systems with 
domain walls, and the conventional z = 2 and 3 for sys- 
tems with conserved and nonconserved order parameter 
respectively have been recovered ^-^]. Here we gener- 
alize to systems with 'intermediate' dynamics (to be de- 
fined below), and/or long-range interactions, and extend 
the method to systems with 0(2) symmetry In this way 
we can confront most of the predictions of the energy 
scaling method. Since the present approach requires ex- 
tended defects, however, it cannot treat the 0(2) system 
in d = 2 spatial dimensions, or the 0(3) system in d = 3, 
for which the relevant topological defects (vortices and 
monopoles respectively) are points. We therefore restrict 
our attention to scalar systems in d = 2, 3 and 0(2) sys- 
tems in d = 3. For these cases, all our results agree with 
the energy scaling predictions, although some subtleties 
arise in the case of the 0(2) system due to infrared sin- 
gularities associated with the far field of the vortex line. 

The paper is organized as follows. Section II deals 
with nonconserved dynamics ('model A'). The method is 
introduced in the simple context of domain walls. The 
extension to vortex lines is then discussed, followed by 
the generalization to long-range interactions. Section 



I. INTRODUCTION 

Coarsening is a characteristic property of a system 
quenched into an ordered phase from a disordered phase, 
and describes the establishment of long-range order over 
ever larger length scales as time increases M. We con- 
centrate here on systems which exhibit the property of 
'dynamical scaling', in which the pair correlation func- 
tion C(r, t) = (<fi(x, t) 0(x + r, £)) of the order-parameter 
field </>(r, t) has the scaling form C(r,t) = f[r/L(t)\. The 
characteristic length scale L(t) generally has a power- 
law dependence on time, L(t) ~ t 1 ^ 2 , sometimes with 
logarithmic corrections, where z is a kind of dynamic ex- 
ponent for the coarsening process. 

Much effort in recent years has been devoted to the 
determination of the exponent z (and any additional log- 
arithmic factors) for a range of different models. Al- 
though the experimentally most relevant case of a scalar 
order parameter has been understood for some time, 
through intuitive arguments based on the scaling assump- 
tion, and through exactly solved limits |J, a general 
approach has been lacking until recently. The 'energy- 
scaling' approach of Bray and Rutenberg filled 
this gap by providing a very general approach based 
on the role of the topological defects nucleated during 
the quench. These defects are simple domain walls for 
a scalar field, vortices/vortex lines for a two-component 
vector field, and monopoles for three-component vector 
field (where we are assuming rotational invariance for the 
last two models, i.e. these models possess 0(2) and 0(3) 
symmetry respectively). The energy-scaling method is 
suitable for any system for which the scaling hypothe- 
sis holds and the dynamics are purely dissipative. These 
include the standard models 'A' and 'B' of the Hohcn- 
berg and Halperin classification jg], describing the sim- 
plest systems with non-conserved and conserved order 
parameter respectively. In the context of phase-ordering 
kinetics, the models arc described by the time-dependent 
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Ill deals with conserved ('model B') and 'partially con- 
served' dynamics, for domain walls and vortex lines with 
both and short-range interactions. Detailed comparisons 
with the predictions of the energy scaling method are 
made at each stage. In particular, the energy scaling pre- 
diction L(t) ~ (tint) 1 / 4 for the conserved 0(2) model in 
d = 3 (or, more generally, d > 3) is recovered. The paper 
concludes with a discussion and summary. 

II. NONCONSERVED DYNAMICS 

A. Domain walls 

The simplest exemplar of the defect relaxation ap- 
proach is the dynamics of an interface (domain wall) 
separating two equilibrium phases. The time-dependent 
Ginzburg Landau (TDGL) equation is a continuum dy- 
namical model for such a process: 

d t( f> = v 2 0-v'(0) , (i) 

where 0(x, t) is the scalar order-parameter field and 
V{4>) is a symmetric double-well potential, e.g. V(4>) — 
(1 — 2 ) 2 /4, whose minima = ±1 represent the two 
bulk phases. 

Consider first a single planar interface separating the 
two phases. The normal to the interface define the z- 
direction. The order parameter depends only on z, and 
satisfies the time-independent version of (|l|), d 2 (j>o/dz 2 = 
V'((j>o), with boundary conditions o (±oo) = ±1. We 
now impose a small periodic perturbation to the inter- 
face in the x-direction with wavevector k 

4>(x, z, t) = 4>o{z) + A<t>i{z) exp(ikx — u>kt) , (2) 

where the amplitude A is small. The relaxation rate u>k 
defines the timescale, = l/u>k> for the relaxation of a 
perturbation with characteristic length-scale 2n/k. Sub- 
stituting this form into (Q), and linearizing in A, gives the 
eigenvalue equation H(f>i = ujk<t>\, with 'Hamiltonian' 

H=-^ + e + V"{M . (3) 

This is conveniently thought of as a quantum mechanical 
Hamiltonian operator. Note that V"(<po) is positive at 
z = ±oo (where 0o = ±1), and negative at z — (where 
0o = 0). For the specific case V(4>) — (1 — 2 ) 2 /4 one has 
V"((j)o) — 30q(z) — 1, which equals — 1 at z = and tends 
to +2 for z — > ±oo. It follows that V"(0o) represents a 
potential well which must have at least one bound state. 
In fact, since k = corresponds to a uniform translation 
of the interface, we know that u>k=o = 0, with eigenfunc- 
tion tpi = d(f>o/dz. Also, since this function has no nodes, 
it must be the ground state. Since the fc-dependence of H 
is simply the additive constant k 2 in this simple model, 
it follows that dcf>o/dz is the ground state for all k, with 
eigenvalue 



uJu = k 2 . (4) 

Higher eigenvalues are separated by a gap from the 
ground state, so any component of the corresponding 
eigenf unctions, in the initial displacement of the inter- 
face, relaxes quickly to zero. For this reason, we concen- 
trate here only on the groundstate eigenfunction, which 
is the 'slow mode'. Neglecting any contribution from 
the other eigenfunctions, @ can be written, correct 
to leading order in the amplitude A, as <p(z, x, t) = 
<Pq(z + Aex.p[ikx — Wfci]), corresponding simply to a mod- 
ulated interface located at z(x,t) = — Aexp[ikx — Wfci]. 
Note that this simple interpretation is a consequence of 
4> = d<po/dz, which will not be true in general. 

Equation (Ji|) has the form u>k = k z with z = 2, from 
which infer the coarsening growth law L(t) - = t 1 / 2 
for this, the simplest of our dynamical models. 

B. Vortex Lines 

If the order-parameter field is the two-component vec- 
tor field c/>(x.,t) of the 0{2) model, the topological defects 
are vortex lines. The analogue of (|l|) for this system is 

8t$= V 2 $-dV/d$ , (5) 

where the 'wine-bottle' potential V((j>) is a function of 
\4>\ only. For definiteness, we will use the form V((j)) = 
(1 — |0| 2 ) 2 /4, with minima on the ground-state manifold 
\cf>\ = 1, although the conclusions arc completely general. 
For this potential @) becomes 

d t $= V 2 $+(l-\$\ 2 )$ . (6) 

This equation has a time- independent solution </>o( r )j cor- 
responding to a vortex line in the z direction, where we 
have introduced the coordinate system x = (r, z) of coor- 
dinates normal and parallel to the vortex. This is a solu- 
tion of (||) satisfying the boundary conditions ^o(O) = 0, 
0o ( r ) — ► r for |r| — > oo, where f is a unit vector. 

By analogy with (||) we now perturb the vortex line by 
adding a periodic modulation with wavevector k in the 
form 

$(r, z,t) = 0o (r) + A0i(r)exp(ifcz - uj k t) . (7) 

Substituting in (^|), and linearizing in A, gives the eigen- 
value equation 

(V 2 - k 2 )$ x + (1 - \$ \ 2 ) $ x - 2(0! • O ) 00 = ~Uk$l , 

(8) 

where V 2 is the Laplacian operator in the r-plane. From 
the same physical considerations that we used for do- 
main walls, we expect a null eigenfunction at k — 0, cor- 
responding to a uniform displacement of the vortex line 
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transverse to its length. Since 4>o(r+a) = <?!>o+a- V r c/>o, to 
first order in a, wc identify a family of null cigenfunctions 
= a • V r </>o, parametrized by the direction of a. For 
convenience we may choose a basis set of two such eigen- 
functions corresponding to a being a unit vector along 
the x and y axes (where r = (x,y)) respectively. The 
two basis functions are then d x (j)Q and d y <po. They are 
orthogonal by symmetry. Substitution into (^) confirms 
that these are null cigenfunctions, i.e. ^k=o = 0. Since 
the fc 2 term in (||) can simply be absorbed into u>k, it 
follows that, just as for domain walls, the smallest eigen- 
value for any fc is exactly given by 



uj k = k 



(9) 



We conclude that the dynamical exponent is again z — 2. 
This agrees with the energy-scaling result. 

For future reference we note that the eigenfunctions 
d x <pQ and d y 4>$ are not normalizable in an infinite two- 
dimensional space, since the normalization integrals are 
logarithmically divergent (see below) . When the normal- 
ization is important, we will take the vortex line to lie 
along the symmetry axis of an infinitely long cylinder of 
radius R. 



C. Long-Range Interactions: Domain Walls 

The above approach may readily be generalized to the 
case where the underlying interactions are long-ranged in 
space. Consider, for example, a ferromagnetic model, in 
which the exchange interaction J(|r — r'|) falls off with 
distance as r — r'|~( d+cr ) for large |r — r'|. The long-range 
character of these interactions is typically only 'relevant' 
(i.e. affects the scaling behavior) if a is smaller than a 
critical value a c (9). In the present context, we will find 
that a c < 2, so we will consider only the case a < 2 here. 
For this case, the term in the Ginzburg-Landau free- 
energy functional generated by long-range interactions 
has the Fourier-space form Flr — (g/2) ^ k fc^kcA-k- 
This can be formally represented in the TDGL dynamics, 
dt<p = —SF/8(f), by replacing the usual Laplacian opera- 
tor by —(—V 2 ) '/ 2 , to give 



d t 6 = V 2 



(10) 



instead of ([j]) . Expanding around a flat domain wall us- 
ing (||), where <fio is now a stationary solution of the full 
equation (fL(J), gives the eigenvalue equation 



(d 2 z -k 2 - g[k 2 - d 2 ]"' 2 - V" 



= -cjk4>i , (11) 



where d z = d/dz. 

The function cj>i = d<po / dz is again a null eigenfunction 
for k — 0, since k = corresponds to a simple translation 
of the interface. The term in k 2 in ( fj"l| ) simply shifts the 
eigenvalue by fc 2 , as before. However, the long-range part 
of (0) modifies the eigenfunction. We can, nevertheless, 



compute the desired small-fc behavior of the ground-state 
eigenvalue perturbatively, using the unperturbed ground 
state eigenfunction d<fio/dz = <fi' {z): 



u> k = k 2 



n, dz <t>> (z) ([fc 2 - diY' 2 - i-4r /2 ) #>(*) 



S^dzW {z)f 



(12) 



The integral in the numerator is conveniently evaluated 
in Fourier space. The function <j)' (z) is sharply peaked 
at the interface, with a width of order £, the interfacial 
thickness, and a peak height of order l/£. Its Fourier 
transform, therefore, is very broad (with width l/£) and 
for q£ <C 1 is equal to 6q(oo) — 0o(— 00) = 2. The integral 
in the numerator in (|l2|) becomes, therefore 
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2lcr/2 



for small fc, provided the integral converges. The latter 
condition requires a < 1. In this regime, the impor- 
tant values of q in the integral are of order fc, so our 
replacement of 4>' (z) by 26 (z) is justified for ^ < 1. 
The integral is then easily evaluated to give k 1+a (up to 
constants) for a < 1. 

For a > 1, the replacement of 4>' a by a delta function 
is no longer valid, since the integral in the numerator in 
( |l2| ) would not converge. The Fourier transform of <fi' 
falls off at q£ w 1, however, to converge the integral. In 
this case (a > 1) the characteristic value of q is of order 
l/£, instead of fc, and we can expand the integrand up to 
order fc 2 for fc£ <c 1. The numerator then becomes of or- 
der fc 2 ^ 1_CT for a > 1. A similar line of reasoning for the 
marginal case a = 1 leads to a fc 2 ln(l/fc£) behavior for 
fc^ < 1. The simple fc 2 dependence for a > 1 just renor- 
malizes the amplitude of the leading fc 2 term in (|l^) . The 
integral in the denominator in ( |l2| ) is a constant of order 

ve 

Putting it all together, retaining only the leading 
small-fc behavior, and discarding constant prefactors, 
gives the dispersion relation 



,1+cr 



a < 1, 



fc 2 ln(l/fc£), 
fc 2 , (7>1. 



£7 = 1, 



(13) 



We deduce that the dynamic exponent is z = 1 + a for 
<7 < 1 and z — 2 for a > 1. The marginal case a = 1 
gives (with fc — > 1/L, ujk 1/t as usual) the growth 
law L ~ (ilnt) 1 / 2 . All these results are in accord with 
previous results based on renormalization group || and 
energy scaling arguments. 
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D. Long-Range Interactions: Vortex Lines 

The influence of long-range interactions on the dynam- 
ics of a vortex line may be discussed in a similar way. 
Taking, as before, V{S) = (1 — |(/>| 2 ) 2 /4, one obtains, 
analogous to (||) and (11), 



g[k 2_ V 2 r /2 + l 

-2(0i ■ 4> )4>o : 



Ao| 2 )0i 



(14) 



The long-range part can be treated perturbatively for 
small fc as in the scalar case. Since the perturbation is 
isotropic, the unperturbed eigenfunctions d x (f>o and d y <fio 
are not mixed by the perturbation. It follows that, to 
leading order in perturbation theory, 



uj k = k + g 



jd 2 rd^ 0] ([fc 2 - V 2 ]-/ 2 - l-Vtr^drfoj 



Jd 2 r[d, 



W0j\ 



(15) 



where the result has been written in a rotationally in- 
variant form, i.e. there are implicit summations over the 
indices i and j . 

We again evaluate the integral in the numerator in 
Fourier space. In the small-fc limit only small wavevector 
q, corresponding to large |r|, will be important. For large 
|r|, (j) -> f, giving d^oj = (<% - f l f j )/\r\. The Fourier 
transform of this quantity is {2i:/q)(5ij — qiqj). Inserting 
this in ( [l5| ) the integral in the numerator becomes, up to 
constants, 



fq 

\l\>l/R I 2 



(Ik 2 



-,21(7/2 



|q| CT ) ~ k a ln(fci?) (16) 



for a < 2, this condition ensuring the convergence of the 
integral at large q. The lower cutoff at q = 1/R, where 
we recall that R is the radius of the system in the plane 
normal to the vortex line, is necessary to regulate the 
logarithmic singularity at small q. 



Now consider the integral in the denominator of (15). 
Using the large-r form of 4>o everywhere gives the re- 
sult J d 2 r/r 2 , i.e. a logarithmically divergent integral. 
At small r, this can be cut-off at the vortex core size 
£, where the assumed large-r form no longer holds. The 
large r cut-off is again R. The need for a large-distance 
cut-off is associated with the non-normali zabil ity of the 
null eigenfunctions, as discussed in section |IIB| . The up- 
shot is that the denominator is of order ln(i?/£), and the 
leading small-fc form of the dispersion relation is 



u k ~ k° 



ln(fci?) 



a < 2. 



(17) 



Some discussion of this result is in order. If we take 
the limit R — > oo in (|l7]), we obtain the well-defined limit 
ujk ~ k a for a vortex line in an infinite system. We ar- 
gue, however, that this is not the appropriate limit in 



which to discuss the implications of ( |17| ) for the coars- 
ening dynamics. In a system with many vortex lines, 
a characteristic scale L can be associated with the line 
density pv (= length of vortex line per unit volume) via 
pv — l/L 2 . Since the far field of a given vortex line is 
screened out on this scale, it is L rather than R which is 
the appropriate cut-off in the coarsening system. Replac- 
ing R by L in (|l7|), and making the usual identifications 
k ~ l/L and lo^ ~ 1/t, gives the coarsening growth law 
L ~ (t/ hit) 1 /' 7 . This agrees again with the predictions of 
the energy-scaling approach. Similar arguments concern- 
ing the nature of the momentum cutoffs will be necessary 
for the conserved 0(2) model in section III B . 



III. CONSERVED DYNAMICS 

A. Domain Walls 

The standard continuum model for the time evolution 
of a conserved scalar field is the Cahn-Hilliard equation 



d t = -v 2 [v 2 0-u'(0)] 



(18) 



For present purposes this may be conveniently rewritten 
in the form 



v 2 (f> + v'(^) + (-v 2 y 1 d t cf>^o 



(19) 



Inserting the form (g) and linearizing in A leads to the 
eigenvalue equation 

/oo 
dz' G k (z - z')Mz')} 
-OO 

+M 2 + i/"(0o)]<M^Ho , (20) 

dq exp[iq(z — z')] 

(21) 



where 



G k (z-z') = 



, 2tt fc 2 + q 2 
exp(— k\z — z'\) 



2k 



is the Green's function for the operator (fc 2 — d 2 ). 

As before, ipi = 4>' is a null eigenfunction for fc = 0. 
For small fc the term in the second square bracket in (|2^) , 
which vanishes for fc = 0, can be treated perturbatively 
to give 



CJfe = fc 2 



Coo d *J-oo dz ' G ^ Z ~ Z ') *o(*0 



(22) 



The integral in the numerator defines the surface ten- 
sion, 7 = dz (0q) 2 « The function 4>o( z ) ac ts 
like a smeared delta function of width £ and strength 2. 
For {;( < 1, therefore, we can replace G k (z — z') by its 
small- argument limit l/2fc in the denominator, to give 
the result 2/fc. In this limit, therefore, 
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(23) 



The dynamic exponent is z = 3, in accord with the ex- 
pected Lifshitz-Slyozov scaling L ~ i 1 / 3 . A more careful 
derivation of this result has been given by Shinozaki and 
Oono @. 

The generalization to 'intermediate' dynamics, or 'non- 
integer derivative models' |Llj ], in which the leading —V 2 
in the Cahn-Hilliard equation is replaced by (— V 2 )^/ 2 
(i.e. by |k| M in Fourier space, instead of fc 2 ) is straight- 
forward. One simply has to insert in ( p2p the Green's 
function for (fc 2 — d 2 )^^ 2 , given by 



G fe (z-z') 



(ig exp[ig(z — z 
2^ 



(fc 2 



y/2 



(24) 



For fi = 2 our previous result is obtained, while for fx = 
one gets Gfc(z — z') = <5(z — z'), giving the nonconserved 
result ui k — k 2 . For general fx, (|24|) gives the scaling form 
G k (z — z') = fc 1_M /(fc|z — z'|). From the discussion of the 
case fi = 2 it is clear that we require the result only in the 
limit k\z — z'\ -c 1, i.e. we need the small- argument form 
of the scaling function f(x). Straightforward analysis 
gives /(0) = const, for fx > 1, and f(x) ~ x* 1 for jti < 1, 
with f(x) ~ ln(l/x) for /i = 1. Thus there is a change of 
behavior at fx = 1. For /i > 1, G/c(z — z') ~ fc 1_M gives 
(with 7 w l/£) w fe - ^~ 1 fc 1+ ' J . For fx < 1, the double 
integral J dz J dz' \z — z'\^~ 1 cj)' (z)(l)o(z') is of order C^ -1 , 
giving the dispersion relation ~ £ _Al fc 2 , while for fx = 1 
one obtains a;/. ~ ^^fc 2 / ln(l/fc£). For convenience we 
summarize these results below, dropping prefactors in- 
volving £: 



fc 1+ ^, 
fc 2 



fx > 1 



ln(l/fcO : 



/Lt < 1 



(25) 



These results have been derived previously by Onuki . 
The dynamic exponent is z = 1 + fx for fx > 1 and z = 2 
for /i < 1. For fx < 1, the exponent is the same as for 
the nonconserved system: the conservation law is 'irrele- 
vant', in accord with a renormalization group argument 

Using u> ~ 1/t and fc ~ l/£ gives the growth laws 
L(t) - fVU+M) f or M > 1, (i/hii) 1 / 2 for /i = 1, and the 
usual nonconserved result L(t) ~ i 1 / 2 for /z < 1. These 
results are again in complete accord with the energy- 
scaling argument. 

Long-range interactions can be included for a con- 
served order parameter in a straightforward way. Omit- 
ting the details, the final result is a rather obvious com- 
bination of (|l2"|) and (52). For general fi it reads 



SZo dz n, tfe' G k {z - z<) d/ {z) 4>> {z>) 



(26) 



where O k is the operator 



O k 



fc 2 + , 9 [(fc 2 - d 2 z y' 2 - i-dly' 2 ] 



(27) 



Note that the numerator and denominator contain the 
information about the interactions and the conservation 
law respectively. It follows that the dynamic exponent 
is z(a,fx) = min(l + a, 2) + max(/i — 1,0). This clean 
separation of the role the interactions and the conserva- 
tion law is mirrored in the 'energy-scaling' approach to 
calculating growth exponents, where the energy depends 
only on the interactions (a) and the energy dissipation 
rate only on the conservation law (/i). 



B. Vortex Lines 

The treatment of vortex lines in a conserved 0(2) model 
follows the same pattern. Consider first the case of 'sim- 
ple' conservation, fi = 2. The starting point is the Cahn- 
Hilliard equation for vector fields: 



(-V 2 )" 1 9^ = V 2 0+(1-|0| 2 



(28) 



Expanding around the stationary vortex solution c/jq us- 
ing (0) gives the eigenvalue equation 



-V 2 



where 



|0o(r)| 2 )] 0i(r)+2(<Mr)-0i(r)) $ Q 



k 2 -uj k \ d 2 r'G k {v~v')$ 1 {v') 



G k {v 



d 2 q exp(iq ■ r) 



(2tt) 2 fc 2 + q 2 



= 0, (29) 



(30) 



is the Green's function for (fc 2 — V 2 ). 

For small fc, the eigenvalue co k can again be calculated 
using first-order perturbation theory. Rotational invari- 
ance ensures, as before, that the unperturbed cigenfunc- 
tions di<poj are not mixed by the perturbation. The result 
can be written in the form 



uj k = fc 2 



Jd 2 r(V$ ) 2 



Jd 2 rfd 2 r' G k (|r - r'|) V0 O (r) • V0 O (r') 
/dWl0o(q)| 2 



/dW/(fc 2 + « 2 )]|<Mq)| 



(31) 



The Fourier-space form is more convenient for present 
purposes. From the result 4>o ~ * r for r 3> £ it follows 
that |^o(q)| 2 -> (27r) 2 /(7 4 for q£ < 1. Using this result 
for all q, and introducing ultraviolet and infrared cutoffs 
l/£ and 1/R respectively as required, the result for ui k 
can be written as 



u k = fc 



jl'iqdq/q 2 



^ R qdq/[q 2 {k 2 + q 2 )] 

,4 MR/0 

* \n(kR) 



(32) 



5 



to leading logarithmic accuracy. 

One again this result requires careful interpretation. 
Taking the limit R — > oo at fixed fc gives u>k — fc 4 , sug- 
gesting the coarsening growth law L(t) ~ t 1 / 4 . The R- 
dependence enters (B2), however, from the requirement 
to cut off infrared divergences associated with the far field 
of the vortex. In the phase-ordering context, the far field 
is cut off at scale L by other vortex lines. This means we 
should replace R by L in (|2|) (see the parallel discussion 
after (|i7|)). With k ~ 1/L and lju ~ 1/i as usual, this 
leads to the coarsening growth law L(t) ~ (tint) 1 / 4 , in 
perfect agreement with the energy-scaling result |||| . 

The extension to intermediate dynamics, controlled by 
an exponent fj,, is again straightforward. We simply re- 
place (fc 2 + q 2 ) by (fc 2 + q 2 y/ 2 in @, with the result 



Kg/0 

ln(fci?) 



(33) 



to leading logarithmic accuracy, for any /i > 0. The 
coarsening growth law becomes L ~ (t lnt) 1 '^**). 

Finally, long-range interactions can be included. The 
expression for u>k combines the numerator from ( |l5| ) with 
the denominator from (|3l|). Replacing (k 2 + q 2 ) by 
(fc 2 + q 2 )^/ 2 in the denominator, and using @ for the 
numerator, gives the result for general fi and a < 2: 



fc <T+/i 



(34) 



Strictly the numerator and denominator generate factors 
ln(afci?) and ln(/3fc_R) (with a 7^ /3 in general) which we 
have cancelled in (p4). The cancellation is strictly valid 
in the limit R — » 00 at fixed fc. In the coarsening con- 
text, where fc ~ 1/L and R ~ L, these factors are of 
order unity so (p4) still holds, to give z = a + \l. 



IV. CONCLUSION 

We have discussed a general method, based on the re- 
laxation of a sinusoidally perturbed topological defect 
(e.g. a domain wall or vortex line), for inferring growth 
laws in phase-ordering systems. The underlying assump- 
tion is that there is a single characteristic length scale 
('dynamical scaling') so that defect relaxation is either 
the sole or dominant coarsening mechanism, or occurs at 
the same rate as the underlying coarsening process. If 
this assumption does not hold, there need be no connec- 
tion between the relaxation spectrum and the coarsening 
exponent. An example where the present approach fails 
has recently been given Jl2[ . 

In this paper we have considered only systems with 
purely dissipative dynamics. These systems can also be 
addressed using energy-scaling arguments, which involve 
equating two independent estimates of the energy dissi- 
pation rate §j4|. In all cases, the results from the two 
different approaches agree. As yet, however, it has not 
proved possible to extend the energy-scaling technique 



beyond purely dissipative systems. The defect relaxation 
method, however, does not suffer from this limitation. 
Indeed, Shinozaki [[l3| has studied the interfacial relax- 
ation spectrum in an incompressible binary fluid, includ- 
ing hydrodynamic effects, and has obtained Uk ~ fc for 
fc — > 0, consistent with the linear growth, L(t) ~ t, of the 
coarsening scale predicted by Siggia Jlij for bicontinuous 
phases. 

Another recent application of this approach is to the 
coarsening of systems which exhibit lamellar structures 
in equilibrium, such as Raylcigh-Benard convective rolls 
as described by the Swift-Hohcnberg equation, or block 
copolymers in the weak segregation regime [jl5j . The dy- 
namics of these systems has attracted much recent atten- 
tion @. 

In the coarsening regime the lamellae do not form 
parallel stripes but rather exhibit a globally isotropic, 
but locally striped, structure with a characteristic length 
scale (as measured, for example, by the typical radius of 
curvature of the stripes) which increases with time, and 
gives the length scale over which the stripes are locally 
roughly parallel. The stripe pattern itself has, in the 
two-dimensional systems studied numerically, topologi- 
cal defects in the form of disclinations. It is not yet clear 
whether the coarsening is described by a single growing 
length scale in these systems, as different measures of this 
scale give different results |l(| . An analysis of the dynam- 
ics of a modulated lamellar structure along the lines pre- 
sented here gives the relaxation rate u>k — (e/256)fc 2 + fc 4 , 
where e is measure of the quench depth |15| . For shallow 
quenches, therefore, one has ~ fc 4 , suggesting a t 1 / 4 
growth at not-too-late times, whereas the growth rate in- 
ferred from the evolution of the structure factor is closer 
to i 1 / 5 |Q. The disclinations, however, have not been 
included in this approach, and it is possible that these 
coarsen more slowly than the interfaces. More work is 
needed to clarify this point. 

In summary, a study of defect relaxation provides a 
simple way to determine growth exponents in coarsen- 
ing systems. This approach should be reliable when the 
coarsening is described by a single characteristic length 
scale. In the cases studied here, the results obtained 
are identical to those derived using the energy-scaling 
method. While the underlying assumptions (dynamical 
scaling) are the same for both methods, the ideas in- 
volved in the present approach are rather simpler, and 
the method is not restricted to purely dissipative dynam- 
ics. The energy-scaling method, on the other hand, is not 
restricted to extended defects, i.e. it can be used for sys- 
tems with point defects or, indeed, no defects at all. 
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